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Abstract 

We derive new duality relations that link the energy of configurations associated with a class 
of soft pair potentials to the corresponding energy of the dual (Fourier-transformed) potential. 
We apply them by showing how information about the classical ground states of short-ranged 
potentials can be used to draw new conclusions about the nature of the ground states of long- 
ranged potentials and vice versa. They also lead to bounds on the T = system energies in 
density intervals of phase coexistence, the identification of a one-dimensional system that exhibits 
an infinite number of "phase transitions," and a conjecture regarding the ground states of purely 
repulsive monotonic potentials. 

PACS numbers: 05.20.-y, 82.35. Jk,82.70.Dd 61.50.Ah 
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While classical ground states are readily produced by slowly freezing liquids in experi- 
ments and computer simulations, our theoretical understanding of them is far from complete. 
Much of the progress to rigorously identify ground states for given interactions has been for 
lattice models, primarily in one dimension The solutions in (i-dimensional Euclidean 
space M'^ for d > 2 are considerably more challenging. Recently, a "collective-coordinate" 
approach has been used to study and ascertain ground states in and for a class of 
interactions . A surprising conclusion of Ref. is that there exist nontrivial disordered 



ground states without any long-range order 



J], in addition to the expected periodic ones. 



Despite these advances, new theoretical tools are required to make further progress. 

Here we derive new duality relations for a class of soft pair potentials that can be applied 
to classical ground states. We consider soft interactions since they are easier to treat theoret- 
ically and possess great importance in soft-matter systems, such as colloids, microemulsions, 
and polymers j^, Q, [t]. The duality relations link the energy of configurations for a pair po- 
tential v{r) to that for the dual (Fourier-transformed) potential. Applications of the duality 
relations lead to some novel results. 

For a configuration = ri, r2, . . . , rjv of ^ 1 particles in volume C with 
stable pairwise interactions [9], f/(r^) = J2i=i j=i'^i^ij) twice the total potential energy 
per particle [plus the "self-energy" 'y(O)], where v{r) is a radial pair potential function and 

= \rj — rj|. A classical ground-state configuration is one that minimizes U{r^). Since we 
allow for disordered ground states, we consider the general ensemble setting. The ensemble 
average of U for a statistically homogeneous and isotropic system in the thermodynamic 
limit is given by 

(f/(r^)) =t;(r = 0)+p / v{r)g2ir)dv, (1) 

where p = limTv^ooy^oo A^/^ is the number density and g2{r) is the pair correlation func- 
tion. It is crucial to introduce the total correlation function h{r) = g2{r) — 1? which decays 
to zero for a disordered system. We consider those stable radial pair potentials v{r) that 
are bounded and absolutely integrable and call such functions admissible. Thus, the corre- 
sponding Fourier transform v{k) in d dimensions at wavenumber k exists, which we also 
take to be admissible, and 

{U{r^)) = v{r = 0) + pv{k = 0) + p v{r)h{r)dr. (2) 
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Lemma. For any ergodic configuration in R*^, the following duality relation holds: 

/ v{r)h{r)dr = j^J v{k)Hk)dk (3) 

If such a configuration is a ground state, then left and right sides of are minimized. 
Proof: Identity ([3]) follows from Parseval's theorem, assuming that h{k) or the structure 
factor S{k) = 1 + ph{k) exists. From ([2]) and ([3]), we see that the both sides of ([3]) are 
minimized for any ground-state structure, although the duality relation applies to general 
structures [lol |. 
Remark: 

Whereas h{r) always characterizes a point pattern, its Fourier transform h{k) is generally 



not the total correlation function o: 
characterizes a Bravais lattice A 



a point pattern in reciprocal space. It is when h{r) 
ll\ that h{k) is the total correlation function of a point 
pattern, namely, the reciprocal Bravais lattice A. 

Theorem 1. If an admissible pair potential v{r) has a Bravais lattice A ground-state 
structure at number density p, then we have the following duality relation for the minimum 

v{r = 0) + ^v{r) = pv{k = 0) + pY^v{k), (4) 
reA keA 

where the prime on the sum denotes the zero vector should be omitted, A denotes the 

n 

reciprocal Bravais lattice [l2j], and v{k) is the dual pair potential, which automatically 
satisfies the stability condition, and therefore is admissible. Moreover, the minimum Umin 
of U for any ground-state structure of the dual potential v{k), is bounded from above by 
the corresponding real-space minimized quantity Umin or, equivalently, the right side of (jlj), 
i.e., 

Umin < Umin = pv{k = 0) + p^v{k). (5) 

keA 

Whenever the reciprocal lattice A at reciprocal lattice density p = p~^{2tt)~'^ is a ground 
state of v{k), the inequality in ([5]) becomes an equality. On the other hand, if an admissible 
dual potential v{k) has a Bravais lattice A at number density p, then 



Umin < Umin = pv{r = 0) + p2_^ v{r), (6) 

reA 

where equality is achieved when the real-space ground state is the lattice A reciprocal to A. 



Proof: The radially averaged total correlation function for a Bravais lattice, which we now 
assume to be a ground-state structure, is given by h{r) = '^n=i ^n^i''^ ~ ^n) — 1, where 
Si{r) is the surface area of a c/- dimensional sphere of radius r, Z„ is the coordination number 
(number of points) at the radial distance r„, and 6{r) is a radial Dirac delta function. 
Substitution of this expression and the corresponding one for h{k) into (I3l) yields v{r = 
0) + X]n=i ^nv{rn) = pv{k = 0) + pYlin=i ^nv{kn), whcrc Z„ is the coordination number in 
the reciprocal lattice at the radial distance Recognizing that J2n=i ^nv{rn) = Yl'reA ^(^) 
(leading to Umin) and J2n=i^'n-v{kn) = t^f fc) yields the duality relation (j4]). However, 

there may be non-Bravais lattice structures Hl] that have lower energy than the reciprocal 



lattice so that Umin < Umin [13] • Inequality ([6]) follows in the same manner as (15]) when the 

ground state of v{k) is known to be a Bravais lattice. 

Remarks: 

1. Whenever equality in relation ([5]) is achieved, then a ground state structure of the dual 
potential v{k = r) evaluated at the real-space variable r is the Bravais lattice A at density 

2. The zero- vector contributions on both sides of the duality relation (jlj) are crucial in order 
to establish a relationship between the real- and reciprocal-space "lattice" sums indicated 
therein jl4|. 

3. We identify below specific instances in which the strict inequalities in ([5]) and ([6]) apply, 
including a theorem and a one-dimensional system with unusual properties. 

Theorem 2. Suppose that for admissible potentials there exists a range of densities over 
which the ground states are side by side coexistence of two distinct structures whose parent- 
age are two different Bravais lattices, then the strict inequalities in ([5]) and ([6]) apply at any 
density in this density-coexistence interval. 

Proof: This follows immediately from the Maxwell double-tangent construction in the U- 
p~^ plane, which ensures that U in the coexistence region at density p is lower than either 
of the two Bravais lattices. 

As we will see, the duality relations of Theorem 1 will enable one to use information 
about ground states of short-ranged potentials to draw new conclusions about the nature of 
the ground states of long-ranged potentials and vice versa. Moreover, inequalities ([5]) and 
([6]) provide a computational tool to estimate ground-state energies or eliminate candidate 
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ground-state structures as obtained from annealing simulations. We will now examine the 
ground states of several classes of admissible functions. 

Admissible functions with compact support.- Recently, the ground states of a class of 
oscillating real-space potentials v{r) as defined by the family of Fourier transforms with 
compact support such that v{k) is positive for < k < K and zero otherwise have been 
studied {2!, Q]. Clearly, v{k) is admissible. Siito {s^ showed that in the corresponding real- 
space potential f (r), which oscillates about zero, has the body-centered cubic (bcc) lattice as 
its unique ground state at the real-space density p = 1/ (8-\/27r^) (with K = 1). Moreover, he 
showed that for densities greater than l/(8\/27r^), the ground states are degenerate such that 
the face-centered cubic (fee), simple hexagonal (sh), and simple cubic (sc) lattices are ground 
states at and above the respective densities l/(6v^7r'^), \/3 / {1Q\/2tc^) , and 1/{8\/2tc^). 

The long-range behavior of the real-space oscillating potential v{r) might be regarded to 
be unrealistic by some. However, since all of the aforementioned ground states are Bravais 
lattices, the duality relation (|4]) can be applied here to infer the ground states of real-space 
potentials with compact support. Specifically, application of the duality theorem in and 
Siito's results enables us to conclude that for the real-space potential v{r) that is positive 
for < r < D and zero otherwise, the fee lattice (dual of the bcc lattice) is a ground 
state at the density \/2 and the ground states are degenerate such that the bcc, sh and sc 
lattices are ground states at and below the respective densities (3\/3)/4, 2/ ^/S, and 1 (taking 
D = 1). Specific examples of such real-space potentials, for which the ground states are not 
rigorously known, include the "square-mound" potential [f (r) = e > for < r < 1 
and zero otherwise] and the "overlap" potential [Sj, equal to the intersection volume of two 
d-dimensional spheres of diameter D whose centers are separated by a distance r divided by 
the volume of a sphere, and thus has support in the interval [0, -D) 15j. The d-dimensional 
Fourier transforms of the square mound and overlap potentials are e2'^^'^ Jci/2{k) / {kjcY^"^ and 
2'^7r'^/^r(l + d/2)J'^i2{k/2)/k'^, respectively, with D = 1. Figure [T] shows the real-space and 
dual potentials for these examples in R^. The densities at which the aforementioned lattices 
are ground states are easily understood by appealing to either the square-mound or overlap 
potential. The fee lattice is a ground state at the density \/2 because at this value, where 
the nearest-neighbor (NN) distance is unity, and lower densities, the energy is zero. At a 
slightly higher density, each of the 12 nearest neighbors contribiutes an amount e to the 
lattice energy. At densities lower than y^, there are an uncountably infinite number of 
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degenerate ground states. This includes the bcc, sh and sc lattices, which join in as ground 
states at and below the respective densities (3a/3)/4, and 1 because those are the 

threshold values at which these structures have lattice energies that change discontinuously 
from some positive value (determined by nearest neighbors only) to zero. Moreover, any 
structure, periodic or not, in which the NN distance is greater than unity is a ground state. 

However, at densities corresponding to NN distances that are less than unity, determi- 
nation of the possible ground-state structures is considerably more difficult. For example, 
it has been argued in Ref. [g] (with good reason) that real-space potentials whose Fourier 
transforms oscillate about zero will exhibit polymorphic crystal phases in which the particles 
that comprise a cluster sit on top of each other. The square-mound potential is a special 
case of this class of potentials and the fact that it is a simple piecewise constant function 
allows for a rigorous analysis of the clustered ground states for densities in which the NN 
distances are less than the distance at which the discontinuity in v{r) occurs [isl . 
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FIG. 1: Left: The localized square-mound potential [v{r) = e = 1 for < r < 1 and zero otherwise] 
and overlap potential [v{r) = 1 — 3r/2 + for < r < 1 and zero otherwise] in M^. Right: The 
delocalized dual square-mound potential v{k) = 7r'^/^J3/2(A;)/(2A;)^/^ multiplied by ir^/G and dual 
overlap potential v{k) = Qi:"^ J^i2{k/2) /k'^ . 

Nonnegative admissible functions. - Another interesting class of admissible functions are 
those in which both v{r) and v{k) are nonnegative (i.e., purely repulsive) for their entire 
domains. The "overlap" potential discussed above is an example. Here we show that the 
dual v{k) = AsirL'{kD /2) / {kDY of the overlap potential for d = 1 [v{r-) = 1 — r/D for 
Q < r < D and zero otherwise] (see Fig. [2]) exhibits rich behavior. For any density p, it 
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FIG. 2: Left: The "localized" overlap potential in M: v{r) = 1 — r for < r < 1 and zero otherwise. 
Right: The corresponding "delocalized" dual potential: ^(A;) = 4sin^(/c/2)//c^. 

can be shown that the unique Bravais (integer) lattice with spacing is a ground-state 
structure j^, [3]. Moreover, using Theorem 1, we can show that for any p = m, where m 
is a positive integer, the integer lattice at reciprocal density p = (27rm)^^ is a ground-state 
structure for for the dual potential v{k); however, at non-integer density p, ground-state 
structures for v{k) are generally non-Bravais lattices, establishing the strict inequality of 
duality relation ([5]) [l3^. This latter result implies that for v{k), the system undergoes an 
infinite number of "phase transitions" from Bravais to non-Bravais lattices over the entire 
density range. This one-dimensional example is interesting in its own right and further 
details about its ground states will be given elsewhere [3]. 

Another interesting example of nonnegative admissible functions is the Gaussian (core) 
potential v{r) = eexp[— (r/a)^] flS^, which has been used to model interactions in poly- 
mers [5]. The dual potentials are self-similar Gaussian functions for any d. The poten- 
tial function pairs for the case d = 3 with e = 1 and a = 1 are v{r) = exp(— r^) and 
v{k) = TT^/^ exp(— A;^/4). Simulations [jj] indicate that at sufficiently low densities in M^, 
the fee lattices are the ground state structures for v{r). For the range < p < 7r-3/2^ fee 
is favored over bcc [19^. If equality in is achieved for this density range, the duality 
theorem would imply that the bcc lattices in the ran ge ( 47r)~^/^ < p < oo are the ground 



state structures for the dual potential. Previous work [181] has verified that this is the case, 
except in a narrow density interval of fcc-bcc coexistence 0.17941 < p < 0.17977 around 
p = TT^^/^ ^ 0.17959. In the coexistence interval, however. Theorem 2 states the strict 
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inequalities in ([5]) and ([6]) must apply. Importantly, the ground states here are not only 



non-Bravais lattices, they are not even periodic 



161]. In M^, the triangular lattices appar- 



ently are the ground states for the Gaussian potential at all densities (but there is no proof) 



and therefore would not exhibit a phase transition. Proposition 9.6 of Ref. [20| enables us 
to conclude that the integer lattices are the ground states of the Gaussian potential for all 
densities in M. 

Completely monotonia (CM) admissible functions. - A radial function /(r) is completely 
monotonic if it possesses derivatives /("^(r) for all n > and if (-!)"/(") (r) > 0. Not 
all CM functions are admissible (e.g., the power-law potential in M"^ is inadmissible) . 
Examples of admissible ones in M."^ include exp(— ar) for a > and l/(r + a)^ for a > 0, 
(3>d. 

Remarkably, the ground states of the pure exponential potential have not been studied. 
Here we apply the duality relations to the real-space potential v{r) = exp(— r) in M.'^ and 
its dual i!{k) = c{d)/{l + P)('^+i)/2 [where c{d) = 2'^7T^'^-^'^/^r{{d + l)/2)], which has a 
slow power-law decay of for large k. Note that ^(A;) is a CM admissible function 

in k'^, and both v{r) and v{k) are nonnegative admissible functions. We have evaluated 
lattice sums for the exponential potential for a variety of Bravais and periodic structures 
in and M.^. In M^, we found that the triangular lattices are favored for all p (as in the 
Gaussian case). If equality in ([5]) is achieved, the triangular lattices are also the ground 
states for the slowly decaying dual potential v{k) = 2tt/(1 + A;^)^/^ for all p. In M^, the 
fee lattices are favored at low densities (0 < p < 0.017470) and bcc lattices are favored at 
high densities (0.017470 < p < oo). The Maxwell double-tangent construction reveals that 
there is a very narrow density interval 0.017469 < p < 0.017471 of fcc-bcc coexistence. The 
exponential potential appears to behave qualitatively like the Gaussian. If equality in ([5]) 
applies outside the coexistence interval. Theorem 1 would predict that the ground states of 
the dual potential v{k) = Svr/ (1 + k^y are the fee lattices for < p < 0.230750 and the bcc 
lattice for 0.230777 < p < oo |22|. 

Conjecture. The Gaussian potential, exponential potential, the dual of the exponential 
potential, and any other admissible potential function that is completely monotonic (CM) 
in distance or squared distance share the same ground-state structures in M'^ for 2 < d < 8 
and d = 24, albeit not at the same densities. For any such potential function, the ground 
states are the Bravais lattices corresponding to the densest known sphere packings [21] for 
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< p < Pi and the reciprocal Bravais lattices for p2 < p < oo, where pi and p2 are 
the density limits of phase coexistence of the low- and high-density phases, respectively. 
Whenever, the Bravais and reciprocal lattices are self-dual {d = 2,4,8 and 24) pi = p2, 
otherwise p2 > pi (which occurs for d = 3,5,6 and 7). 



This conjecture is bolstered by the recent work of Cohn and Kumar 20], who rigorously 



proved that certain configurations of points interacting with CM potentials on the surface 



23| 



of the unit sphere in arbitrary dimension were energy-minimizing 

In summary, we have derived and applied duality relations to help quantify and identify 
ground states for pair potentials that arise in soft-matter systems. Elsewhere, we will apply 
the duality relations to a broader category of functions beyond the pure Gaussians that are 
self-similar under Fourier transformation and will show that our formalism be extended to 
obtain corresponding duality relations for potential functions that also include three-body 



and higher-order interactions ISj. 
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